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Abstract 

In this study, we examine decoherence of a spin qubits system coupled independently to counted 
spin chain with a 1/r 2 interaction by using influence functional. We also examine the time evolution 
of density matrix numerically when environment is Gaussian noise. 
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I. INTRODUCTION 



Among various proposals for quantum computations, quantum bits (qubits) in solid ma- 
terials, such as superconducting Josephson junctions^, and quantum dots 2 - 1 ^, have the ad- 
vantage of scalability. Such coherent two level systems constitute qubits and a quantum 
computation can be carried out as a unitary operation applied to many qubit systems. It 
is essential that this quantum coherence be maintained during computation. However, de- 
phasing is hard to avoid, due to the interaction between qubit system the environment. 
The decay of off-diagonal elements of the qubit density matrix signals occurrence of dephas- 
ing. Various environments can cause dephasing. In solid systems, the effect of phonons 
is important^. The effect of electromagnetic fluctuation has been extensively studied for 
Josephson junction charge qubits 6 . 

For spin qubit system, the fluctuations of the nuclear spins of impurities can also be a 
cause of dephasing. It has recently been shown experimentally that the coupling between the 
spin of an electron in a quantum dot and the environment is very weak"*^. Decoherence of 
nuclear spins is experimentally examined, in which the source of decoherence is dipole-dipole 
interaction.— For this reason, the dephasing time of a spin qubit conjectured to be very long. 
However, both donor impurities and nuclear spins in semiconductors 8 - have been suggested 
as possible building blocks for feasible quantum dots architectures. The proposal based on 
experimental findings of quantum computation by using an Si 29 array is another possibility 8 . 
For a spin qubit system in 2D, the coupling with neighbor spin chain can cause a dephasing. 
When an interaction is between the qubits themselves, it can in principle be incorporated 
into the quantum computer Hamiltonian, although this would lead to more complicated 
gate sequences. Therefore it is instructive to analyze the error introduced by ignoring some 
of these interactions, as done in the case of dipolar coupled spin qubits^ and spin-orbit 
interaction.— For a quantum spin chain with a 1/r 2 interaction, an exact expression of 
the dynamical correlation function has been obtained^. Using this expression, we consider 
the case in which the qubit system is coupled to spin of the spin chain. We examine the 
relaxation phenomena of a spin qubits array which is coupled to a spin chain with long-range 
interactions. Integrating over the spin chain variables, we obtain the influence functional 
of the qubit system.—^ Using this influence functional, we examine the dephasing of the 
qubits density matrix. In the present study, we especially concentrated on the effect of spin 
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flip process. We examine the zero-dimensional qubit and one- dimensional qubits system. 
The spin flips process lead to oscillation self-excitation.— 

II. HAMILTONIAN 

We examine the Hamiltonian of counted spin chain and qubits system. The type of pin 
chain and qubit interaction is XXZ, (A±, A zz ^ 0). 

H spin = J^(-ir- m [d(n-m)] 2 S n -S m , (1) 

n,m 

H qb = Y,^oI?, (2) 

i 

H int = ^Nh 2 [~A ± (S?I? + S?I?)+A zz S*m. 

i 

(3) 

This Hamiltonian H sp i n represent inverse square interacting spin chain system, and is counted 
Haldane-Shastry model. Here, d(n)=(N/ir)sm(im/N), Sj,(j = x,y,z) is spin operator at 
site i of spin chain, if ,(j = x, y, z) is qubit operator at site i, J is strength of interaction of 
spin chain system, 7 at is geomagnetic ratio of qubit. Above model, the spin chain and zero- 
dimensional qubit interacts with contact interaction. We also examine the one-dimensional 
qubit system with contact interaction, for this case the effect of indirect interaction appears. 
And, we scale the length by the lattice length a. 

III. SINGLE QUBIT SYSTEM 

First, we consider a single qubit and one- dimensional spin chain system. Thus, we ex- 
amine the effect of direct interaction by influence functional method.— 1 ^ The interaction 
Hamiltonian is as below, 

H^t = i N h 2 [\A ± (S*r + SIP) + A ZZ S Z I Z ]. (4) 

We integrate about the qubit system except the 0-th site spin. This lead density matrix 
of spin system, 

p(I f z+ ,lL) = / [dl z+ ][dl z .}exp(-(l qb [l z .] - I qb [I z+ ]))F[I z+ ,I z „], (5) 
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the influence functional is 



F[I Z+ ,I Z _] = J [dS +i ][dS. 
(J(S +i (t)-S_ i (t))p(S i (0) J S i (0)) 



exp{|(J[S + ]-/[S_])}, (6) 

where, = f^hAI z . For following discussion, we choose the action of system, I[S] = 

io[S] + 4(^2, S^ ], thus the unperturbed is given by, 

\2 ft rtl 



J [S] = ^^jf j 1 dt 1 dt 2 Al z St{t l )A^\t 1 ,iM,i)Sm 



Al(St(t 1 )A^ p -\t u z,t 2 ,t)S-(t 2 ) + S7( tl )A 00 + p '\t ui ,t 2 ,z)Snt2)) 



+ t— dx dtn ■ (— x — ) 7 
47T7 7 v &r <9r 

Where # = 27m and Aoop{ti,i,t 2 ,j) is the free propagator of environmental system at zero 
temperature which is defined on closed time path and has four components, 

Aoop(i,ti,j,t2) 

A oo + (*> *i> i> fe) A oo~ (h h, j, t 2 y 

(8) 

For the Berry phase term, because n is odd, the environment is a one- dimensional spin-1/2 
system.— In the present model, the spin chain is a gapless solvable model (a special point in 
spin chain systems). Therefore, we can integrate over the spin chain degree of freedom and 
the Berry phase term does not make system to be distempered. Introducing the incoming 
interaction picture for the environment system we can easily verify that Eq. (JBJ) turns out,— 

F[I Z+ ,I Z _] = expH^^ f f 1 dhdt 2 

4 Jo Jo 
Al z {I z+ {t 1 )A+ + (i,t 1 ,i,t 2 )I z+ (t 2 ) 

+ J 2 _(ii)A 00 ~(2,ti,z,£ 2 )/ 2 _(t 2 ) 

- J 2+ (ti)A+"(i,ti,2,t 2 )4_(t 2 ) 

- I z _(t 1 )A 00 + (t,t u t,t 2 )I z+ (t 2 )) 
-Alm+ih) A++(i, t 1} i, t 2 )I++(t 2 ) 



+ /+_(ti) A^j-(i, ti, i, t 2 )/+_(* 2 ) 

- / ++ (* 1 )A+-(i,t 1 ,i,t 2 )/+_(* 2 ) 

- 7 + _(i 1 )A - + (i,t 1 ,i,t 2 )7 ++ (t 2 )) 

(/_+(*!) A++(*, *!, i, * 2 )/_+(t 2 ) 
+ I—(t 1 )A--(i,t 1 ,i,t 2 )I__(t 2 ) 

- I_ + (t 1 )A+-(i,t 1 ,i,t 2 )I^(t 2 ) 

- I_-(t 1 )A£-(i,t 1 ,i,h)I- + (t 2 ))]. 

(9) 

For the convenience we change of the coordinate, 

Vz = {I z+ + h-)l%i z = (I z+ - I z -)/2. (10) 

The rj and £ are called sojourn and blip. In terms of this variable, density matrix is described 
by as follows 

p(r] z = 1) = | t><T \,p(Vz = -i) = | !><! |, 

p(6 = i) = IT><Hp(& = -i) = U><U 

Another variable are defined by 

t? + = (J ++ + 7 + _)/2, e + = (/++ - /+-)/2, (11) 

77_ = (7_+ + 7__)/2, = (7_ + - 7__)/2. (12) 
In terms of above variables the elements of density matrix are expressed as 





= 1) = 


1 1><! \,p(v+ = 


-i) = 


1 IX 1 


p(£+ 


= 1) = 


ir><Ti,p(e+ = 


-i) = 


1 1><I !■ 


p(v- 


= 1) = 


| |><T = 


-i) = 


1 XI 1 




= 1) = 


1 1><I I,p(£- = 


-i) = 


1 T><T 1- 



Therefore, the below equations are hold for these new variables, 

Vz = t+ = -£-,tz = v+ = -v- ( 13 ) 
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In this case, the influence functional is expressed as 

\2 ft r ti 



F[r], £] = expH^- J J'dhdh 

A 2 z Mti)G R (h,i,t2,i)v(t2) 
+ r ? (t 1 )G ,A (t 1 ,^,t 2 ,2)e(t 2 ) 



- ati)G K {h,i,t2,mt2)} 

+ £(fi)G A (ti,M 2 ,i)»7(f2) 

- r ) (t 1 )G K (t 1 ,i,t 2 ,i)r ] (t 2 )}) 

+ £(ti)G A (*i,M2,iM* 2 ) 

- riitjGK^i^iWh)})} 

= exp[—i —^^- I I dtidt 2 
4 Jo Jo 

{^(AlG^ i, t 2 , i) + A\G A {ti, i, t 2 , i))rj(t 2 ) 
+ 77(*i) (A 2 zz G A (t u i, t 2 , i) + A 2 ± G R (t u i, t 2 , i))£(t 2 ) 

- 4U(*i)g*(M,M)£(* 2 ) 

- 2A 2 zz r ) (t 1 )G K (h,i,t 2 ,i)r ] (t 2 )}}, 

(14) 

where G R (ti, i, t 2 , j), G A (ti, i, t 2 , j) and G K (ti,i,t 2 ,j) are retarded Green's function, ad- 
vanced Green's function and Keldysh Green's function. For above integral equation, after 
we slice the time and take a difference, we get the differential equation for density matrix is 
given by, 



dph(t) = ~[H qb ,pm-^Nh) 2 £d tl 

\- /.;./,). a 2 G A a.t.i.u\ + \- c'Ui.i.i.i^ \ 

Pb(h) 



-A 2 zz G K (i, t, i, h), A 2 zz G A (i, t, i, h) + A 2 ± G R (i, t, i, h) 

A 2 zz G R (i, t, i, h) + A\G A (i, t, i, h), -2A 2 ± G K (i, t, i, h) 



(15) 
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where [A, B] is AB - BA and p h (t) is 

Vz {t) = 1, tfe(t) = -1 
£,(*) = 1, £.(*) = -1 
Next we choose the representation of density matrix for the spin diagonal case, 



Pb(t) 



/ 



± zz' Ji ± 



A lz+ A 



G K (i,t,i,t x ) 
HG R (t,t,t,t 1 ) + G A (t,t,t,t 1 )) 



\ 



«jg^ (G A (i,t,i,t 1 ) - G R {i,t,iM)) 



\ 



42 42 „ 



. where the density matrix, p s (t) is represented as 

/|T(*)><T (*)I + U(«)><i 
IT (*) ><I (*)|, 

\i(t) ><T 
VlT(t) ><T (t)\-\i(t) ><i(t)\J 



Ps{t) 



(16) 



(17) 



(18) 



The trace of density matrix decreases with time. Another diagonal element show different 
behavior. For spin flip process, another diagonal element increases with time, this represents 
self-excitation. The one of off-diagonal element become decoherence. Another off-diagonal 
element shows oscillation where modulation of signal occurs. 



IV. ONE-DIMENSIONAL QUBITS SYSTEM 

Next, we examine 1-dimensonal qubit system by using influence functional.— In this case, 
the effect of indirect interaction appears. The interaction Hamiltonian is as below, 

H int = iNh 2 Y^A L {ltS: + irS?) + A zz IlSt\. (19) 

i 

We integrate about the qubit system except i-th site spin. This lead density matrix of spin 

system, 

~i zi+ {t)=i f zi+ ,i zi ^(t)=i f zi _ 



[dl zi+ ] [dl zi _] 



W(o)=/ z J l+ ,/ ZI+ (o)=/ 2 ' i+ 



exp(^(/ ?6 [4 i _] - I qb [I zi+ )))F[I zi+ , I zl _], (20) 



the influence functional is 



F[I zi+ ,I zj _] = J[dS +i ][dS^)5{S +i (t)-S^{t)) 



p(S J (0),S,(0))exp{^(/[S + ]-/[S_])}, 



where, I q b[I] = J hAIiZ. For following discussion, we choose the action of system, I[S] = 
Jo[S] + Iint[Izi, S z i°], thus the unperturbed is given by, 

7 [ S ] = ( A ^ hf f r dt l dt 2 Y,smKi;\ti,iM,3)s z 3 {t 2 ) 

4 J J i.j 

+ iA± l Nh) \ st(t 1 )A+ -\t^t 2 ,j)Sj(t^ 

+ ^JdxJ o dtn.(-x-) (21) 

where 9 = 27m and AQ 0p (ti, i, t 2 , j) is the free propagator of environmental system at zero 
temperature which is defined on closed time path and has four components. For the Berry 
phase term, because n is odd, the environment is a one- dimensional spin-1/2 system.— In 
the present model, the spin chain is a gapless solvable model (a special point in spin chain 
systems). Therefore, we can completely integrate over the spin chain degree of freedom and 
the Berry phase term does not make system to be distempered. Introducing the incoming 
interaction picture for the environment system we can easily verify, and the equation turns 
out, 

F[I zl+ ,I zt _] = expH^£ f I* dt x dh 

4 Jo Jo 
A 2 zz {I zi+ {t 1 )A++(i, t t ,j, t 2 )I zj+ (t 2 ) 

+ J Z4 _(ti)A 00 _ (z, t u j, t 2 )Izj-(t 2 ) 

- 7 ri+ (ti)Aj,~(i, ti,j, h)I zj -(t 2 ) 

- J 2l _(ti)A 00 + (i, t x ,j, t 2 )I zj+ {t 2 )) 
^Al(I +i+ (t 1 )A++(i,t l ,j,t 2 )I +j+ (t 2 ) 

+ I +i _(£i)Aoo~ (i, h,j, t 2 )I+j-(t 2 ) 

- I +i+ (ti)A^(i, t 2 )I+j-(t 2 ) 

- J +i _(ti)Ao + (z, ti,j, t 2 )I +j+ (t 2 )) 



^Al(I_ i+ (t 1 )A+ + (i,t 1 J,t 2 )I. j+ (t 2 ) 

- /_ i+ (t 1 )A+-(i,* 1 ,j,* 2 )/_ j _(t 2 ) 

- /_ i _(t 1 )A 00 + (2,t 1 ,j,t 2 )/_, + (t 2 ))]. (22) 
For the convenience we change of the coordinate, 

Vzi = {hi+ + Izi-)/2, U = (/«+ - hi-)/2. (23) 

The r/ and £ are called sojourn and blip. In terms of this variable, density matrix is described 
by as follows 

p(Vzi = 1) = I U><U \,p(Vzi = -1) = | li><U |, 

pfe* = i) = I h><U \,p(U = -i) = I U><h |. 

Another variable are defined by 

^ = (7 +t+ + /-h_)/2, = (7 +l+ - /-h-)/2, (24) 

^ = (7_ i+ + J_i_)/2, = (7_ i+ - 7_,_)/2. (25) 
In terms of above variables the elements of density matrix are expressed as 



p(v+i = 


i) = 


1 u><u 


l,p(»7+i = 


-i) = 


1 ii><Ti 1 


p(£+i = 


i) = 


1 Ti><Ti 


I,p(£h = 


-i) = 


1 li><U 1 


p(v-i = 


i) = 


1 U><h 


= 


-i) = 


1 h><U 1 


p{U = 


i) = 


I |i><|i 


1 » » = 


-i) = 


1 Ti><T< 1 



Therefore, the below equations are hold for these new variables, 

Vzi = £+i = -f-i, izi = V+i = ~V-i- (26) 
In this case, the influence function is expressed as 
= expH^^ f pdt.dt, 

{Uh){A 2 zz G R {t u i, t 2 , j) + A 2 ± G A (h, i, t 2 , j)) Vj (t 2 ) 
+ ^(Al^^iM.j) + A\G R {t u iM,3)Mh) 
- A 2 zz ^i(ti)G K (ti,i, t 2 , j)£,j{h) - 2A 2 zz r ]i (t 1 )G K (t 1 ,i,t 2 ,j)r ]j (t 2 )}}, 

(27) 



where G R (t\, i, t 2 , j), G A (ti,i,t 2 ,j) and G K (ti,i,t 2 , j) are retarded Green's function, ad- 
vanced Green's function and Keldysh Green's function. For above integral equation, after 
slice the time and take difference, we get the differential equation for density matrix is given 
by, 

dpb(i,j,t) ir . , i. 2 /"* 



dt 



A 2 zz G R (i,t,k,t 1 ) + A 2 , G A (i,t,k,t 



—Al z GK(i, t, k, ti) 



o 

A 2 zz G A (i, t, k, h) + A{G H {i, t, k, h 
-A\G K {i,t,kM) 



k 



Pb(h,k,j) 



(28) 



Next we choose the representation of density matrix for the spin diagonal case, 



dp s (i,j,t) 
dt 



/ -^4^G*(i,f,Mi) 

£^{G R (i, t, k, h) + G A (i, t, k, h)) 
iiA "- A ^ (G A (i, t, k, h) - G R (i, t, k, h)) 



\ 



■ ^lG K (t,t,k,t 1 ) 



p B (k,j,ti) 



J 



(29) 



The trace of density matrix decreases with time. Another diagonal element shows different 
behavior. For spin flip process, another diagonal element increases with time, this represents 
self-excitation. The off-diagonal element shows oscillation where modulation of the signal 
occurs. The self-excitation is appears in poor man scale for Kondo effect. The strong 
coupling limit is self-excitation. Above behavior is simple RG flow. 

We examine the pure dephasing event. Because the propagator of qubit has no time 
dependence, the blip state and sojourn state does not change. Therefore when we choose 
the initial condition of qubit density matrix to be coherent state, such as, 



p(t = o) = n±(ITi><UI±Ui><Ti|), 



(30) 



the time evolution occurs only off-diagonal channel. In addition to that, even if we start 
from off-diagonal state, the interaction Hamiltonian does not contain the spin flip process. 
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Then the blip state does not change. So we can take the ^(t) = &(= ±1) and r)i(t) = for 
all t. This situation leads to exact expression of dephasing rate as follows, 



f i zi+ (t)=if i+ ,i zi .(t)=i f zi _ 
P&,t)= / 

Jl zi+ {0)=I f zi+ ,I zi+ (0)=I f zi+ 



exp[iAt - i {Azz ] Nh)2 v fdt f 1 dU^it^M.m 

= exp(iAt) x detit- r" 4 ^ / d *i / ' dt 2 G K (t 1 ,t,t 2 ,j)} 

16 Jo Jo 

= exp^At-Trlogfl- /^ 7 ^ 2 fdt, P dt 2 G K {t u i, t 2 ,i)}}, 

-Lb Jo Jo 



(31) 



where we neglect the effect of direct interaction. Using the analytic expression, the time 
evolution of off diagonal density matrix element is given by 

- In8fy>(& = ±l,i) 

\2 ft r tl 



Trlog[l + {Azz ^ h)2 J dtj' dhiSfalSfa)}] 
£log[l + {AzZ lf f j/ 1 ^ dt 2 {S z ( tl ),S z (t 2 )}] 
Llog[l+ (Azz ^ h)2 dw{S z (w),S* (-w)} 

^ J — oo 



&m(ut/2) 
Uj2 



(32) 



where {Sq, Sq} is symmetries correlation function, defined by {A, B} = AB + BA because 
Keldysh Green's function is symmetries correlation function for the base of present Hilbert 
space. Then, 

A.M = (1 + £*.{*<«,),$!(-»)} (^p) V (33) 

Above result is exact. Thus, for infinite number of qubit system the dephasing rate becomes 
infinity. Next we evaluate the quality factor. The quantum error correction code rate is given 
by AL, where A is single qubit coherence time and L is number of qubit. The decoherence 
rate is given by ln(l + -^-t)j where T 2 is decoherence rate by single qubit. Therefore the 
quality factor is given by Q = At _ 1 . At t — > oo, Q becomes infinity, thus spin quantum 

ln(l+V*) 

computer is scalable. 
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V. NUMERICAL RESULTS 



For Gaussian noise, we obtain the numerical results. The initial conditions are 
I{t),a x ,o y ,G z = 1. The coupling constants are A zz = 1.0, A ± = 0.1. The off-diagonal 
components show oscillation. The trace shows decoherence. The diagonal component in- 
creases. The purity shows oscillation. The definition of purity is given by Tr(p 2 (t)) 

For SU(2) symmetric case, (A zz = A±), the equation for density matrix is given by as 
follows, 



dl(t) 
dt 
da z (t) 
dt 

dt ~ V j 
da z (t) 



1 /"* 

--( lN hA) 2 / dhCit-h)^ 

4 Jo 



dt =°- ^ 

By using above equations, the equation for a y (t) is 

^)=-^ W+! A( 7 ^/ dtlC(t - tl ,^ ,35, 
The numerical results for SU (2) symmetric case is ginven as Fig. 6-10. 



VI. CONCLUSION 



In summary, we had examined decoherence of a spin qubits system coupled with a spin 
chain. The examined decoherence is the case of single qubit system and one-dimensional 
qubits system. We obtained the differential-integral equation for general initial condition. 
The trace of density matrix decreases with time. Another diagonal element shows different 
behavior. For spin flip process, another diagonal element increases with time, this represents 
self-excitation. The off-diagonal element shows oscillation where modulation of the signal 
occurs. Decoherence without trace conservation occurs. At thermodynamic limit, quality 
factor becomes infinity, thus spin quantum computer is scalable. We also examine the nu- 
merical calculation for Gaussian noise. The purity is plotted, this quantity shows oscillating 
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FIG. 1: Real part of I(t) 



behavior. 
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FIG. 3: Real part of cr, 
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FIG. 5: purity 
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FIG. 6: Real part of I(t) 
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FIG. 7: Real part of a, 




17 




FIG. 9: purity 
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